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We discuss on an example a general mechanism of apparition of anomalous scaling in scale invari- 
ant systems via zero modes of a scale invariant operator. We discuss the relevance of such mechanism 
in turbulence, and point out a peculiarity of turbulent flows, due to the existence of both forcing and 
dissipation. Following these considerations, we show that if this mechanism of anomalous scaling is 
operating in turbulence, the structure functions can be constructed by simple symmetry consider- 
ations. We find that the generical scale behavior of structure functions in the inertial range is not 
self-similar S„(£) oc f> n but includes an "exponential self-similar" behavior S n {£) oc expect -1 £ a ] 
where a is a parameter proportional to the inverse of the logarithm of the Reynolds number. The 
solution also follows exact General Scaling and approximate Extended Self-Similarity. 

PACS numbers: ll.30.-j — Symmetry and conservation laws - 47.27.-i Turbulent flows, convection, and 
heat transfer - 47.27.Gs Isotropic turbulence; homogeneous turbulence - 47.27.Jv High-Reynolds-number 
turbulence 



Resume 

A partir d'un exemple, nous discutons un mecanisme 
general de production de lois d'echelle anormales, dans 
un systeme invariant d'echelle. Ce mecanisme repose sur 
l'existence de valeurs propres nulles pour un operateur 
invariant d'echelle. Nous discutons ensuite la pertinence 
de ce mecanisme en turbulence, en soulignant une par- 
ticularite des ecoulements turbulents, liee a la coexis- 
tence d'un forgage et de la dissipation. En utilisant 
ces considerations, nous montrons que si ce mecanisme 
s'applique a la turbulence, alors on peut construire 
les fonctions de structure par de simples arguments de 
symetrie. On trouve que le comportement generique des 
fonctions de structure n'est pas auto-similaire S n {£) oc 
£^ n mais inclut un terme "exponentiel auto-similaire" 
du type S n {£) oc exp[£ n a _1 £ Q ], ou a est un parametre 
inversement proportionnel au logarithme du nombre de 
Reynolds. La solution satisfait egalement rigoureusement 
la propriety de Similarite Gencrale et approximativement 
la propriete d'Autosimilarite Etendue. 



I. INTRODUCTION 

Anomalous scaling usually refers to situations where 
the observed scaling exponents deviate from their nat- 
ural dimensional scaling. As a consequence, it is often 
associated with non-simple scaling behavior of the corre- 
lation functions: the scaling exponent of the 2n th order 
correlation function is not n times the scaling exponent 
of the 2 nd order correlation function. A good example of 
this situation is the scaling of the so-called velocity struc- 
ture functions in 3D homogeneous isotropic turbulence. 
Dimensional considerations based on energy conservation 
PQ predict that the scaling exponent of the n th structure 



function should be £„ = n/3. Observed scaling expo- 
nents substantially deviate from this prediction (see e.g. 
';2j for a review of recent experimental results). Several 
general mechanisms of apparition of anomalous scaling 
have been identified. Among them, two mechanisms de- 
serve special interest: existence of additional essential pa- 
rameters which cannot be eliminated from the problem 
(this is the phenomenon of incomplete self-similarity, re- 
viewed by [3]) and existence of non-trivial zero modes of 
the closed equations satisfied by the correlation functions 

IS- 

Our aim is to study some consequences of these two 
form of anomalous scaling in situations where the sys- 
tem is scale invariant. Specifically, we show on a simple 
example how zero-modes give rise to a self-similar be- 
havior of second kind. We then come back to the case of 
turbulence to show how this feature can be used to de- 
termine general properties of the scaling exponents and 
the scale behavior of the structure functions. The main 
conclusion of the present study is that the generic scale 
behavior of structure functions in the inertial range is 
not self-similar S n {£) oc £"* n but "exponential self-similar" 
S„(£) oc exp[^ n a _1 £ Q ] where a is a parameter propor- 
tional to the inverse of the logarithm of the Reynolds 
number. The self-similar shape is recovered only in the 
limit of infinite size inertial range (infinite Reynolds num- 
ber). 



II. ANOMALOUS SCALING IN SCALE 
INVARIANT SYSTEMS 

The goal of the present Section is not to prove rigorous 
results, but to illustrate a possible mechanism of genera- 
tion of anomalous scaling in turbulence, and explore its 
consequences. In particular, the discussion about appli- 
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cation to turbulence should be taken as rather specula- 
tive, since we do not derive any exact equations. The 
main goal of the present Section is to give justifications 
for the symmetry approach presented in the next Section, 
which will provide some new and quantitative results on 
structure functions in turbulence. 



A. An example 



scale smaller than L, stationary solutions can be found 
by using the scaling properties of F (Section 2.1) and 
balancing the l.h.s. and the r.h.s. of 0: 

£„F„(ho t (4) 

which predicts another power-law solution cf> cx £ 1+h ", for 
^< L. 



Consider the differential equation: 

d t ^ + d e( t>-^ = F(£), (1) 

where t is the time, £ is the scale, is a scale-dependent 
physical quantity (e.g. a structure function, or a corre- 
lation function), F an external forcing and C a constant. 
We further assume that the scale range is bounded by 
some larger scale L, at which <fi obeys the boundary con- 
dition cf>(£ = L) = 0, and that, for scales much smaller 
than the larger scale, F follows a power-law F(£) cx £ h " . 
This property can be obtained for example if F is a Gaus- 
sian random function with a power-law spectrum. 

1. Scale invariance 

We note that, in absence of forcing, (F = 0), the equa- 
tion is invariant under the family of transformations 
S h (A): 

S h {\):£^\£, -> \ h 4>, t^-Xt, (2) 

for arbitrary A and h. The physical system described by 
eq. without forcing is then said to be scale invariant. 
Several important properties of scale invariant systems 
are listed in Section 3.1. 



4- Incomplete similarity 

The previous dimensional analysis ignores the influ- 
ence of the largest scale L. In fact, the solution must be 
influenced by the scale L, would it be only to satisfy the 
boundary condition <j)(L) = (which is not satisfied if 
the solution is only a power law!). A correct dimensional 
analysis thus imposes: 

cf>(£)=£F(£)G(£/L), (5) 

where G is a function which cannot be determined by di- 
mensional arguments. A modification of the dimensional 
power law solution 4> cx £ 1+h ° can then be observed if the 
variable L is essential and if G obeys incomplete similar- 
ity in L, i.e.: 

G(£/L) ~ (£/L)- A , L > 1, A > 0. (6) 

Since A is positive, the function G tends to infinity when 
L tends to infinity. In other words, the size of the sys- 
tem is always relevant and cannot be eliminated from the 
problem. In such a case, 4> follows an anomalous scaling: 
(f> cx L A £ 1+h °~ A . The value of A can be found by solving 
exactly the equation 0. 

5. Exact solution of the complete problem 



2. Zero-mode 

In absence of forcing, stationary solutions of eq. 
satisfy the homogeneous equation: 

9^-^ = 0. (3) 

The solution is <fi = A£^ , where A is to be chosen accord- 
ing to the boundary condition, namely: 4>(L) = AL^ = 0, 
i.e. A = 0. In absence of forcing, there are therefore no 
non trivial stationary solutions of eq. . The situation 
is modified when the forcing is taken into account. 

3. Dimensional analysis 

Before solving the forced equation, let us perform a 
simple dimensional analysis of eq. Q. We note that for 



The stationary solution of eq. is: 

r l 

(j) = Af+f / p-^F{p)dp. (7) 
Jo 

The constant A is found by imposing the boundary con- 
dition at £ = L: 

A = - [ p^F(p)dp~ L 1+h °-t. (8) 
Jo 

The leading contribution to (f> for £/L <gC 1 then depends 
on the sign of A = 1 + ho — (. When it is negative, <p 
just obeys the dimensional scaling, and the contribution 
due to L tends to zero as L — > oo. This is the case of 
complete similarity, with no anomalous scaling. When 
A is positive, the dominant contribution is , which 
tends to infinity with L. This is the case of incomplete 
similarity described previously. 
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6. Summary 

This simple example illustrates the general mechanism 
of anomalous scaling discussed in [4| in the case of a fi- 
nite size, scale-invariant system: a zero mode, with a 
generic non-dimensional similarity exponent, and which 
does not satisfy the boundary condition, becomes rele- 
vant in the presence of an external forcing, which also 
imposes a dimensional, normal scaling exponent. The 
dominant contribution to the solution then depends on 
the relative value of the two scaling exponents, and in 
particular, on the behavior of the limit L — > oo. 

7. Implications on scaling exponents 

In the case discussed previously, non-anomalous expo- 
nents are given by the exponents of the zero mode of a 
scale invariant differential operator. It then makes a lot 
of sense to investigate the possible constraints on scaling 
exponents dictated by the scale symmetry. This fully jus- 
tifies computation of scaling exponents based on multi- 
fractal theory [5j, conformal theory [6], or more recently 
on a principle of exponent relativity developed in [7J [5] . 
Such approaches restrict more or less severely the pos- 
sible shape of non-anomalous scaling in scale invariant 
systems. Restrictions are very weak in the case of the 
multi-fractal or conformal theory, while they are rather 
strong in the exponent relativity approach: in that case, 
the possible generic exponents fall into three classes de- 
pending on the properties of the rarest events. Within 
each class, the shape of the scaling exponents as a func- 
tion of the order of the correlation function is entirely 
determined by the scale symmetry, and only depends on 
three parameters. 

Transition from anomalous to non-anomalous scaling 
can also be constrained within the present model. Ob- 
viously, the value of the scaling exponents cannot be 
affected by modifications of external parameters, such 
as the strength of the force. The transition to non- 
anomalous scaling can only be obtained by a modification 
of the shape of the differential operator, i.e. by variation 
of internal parameters such as dimensionality, or isotropy 
properties. In the case of the white advected passive 
scalar, for example, it has been proved [§] that the non- 
anomalous behavior disappears in the limit where the 
(topological) dimension of the system is infinite. This 
kind of asymptotic behavior could well be found in a 
large variety of systems, such as turbulence. 

B. The case of turbulence 

The Navier-Stokes equations are: 

d t u + (u ■ V)w + VP = iA7 2 u + F, V • u = 0, (9) 

where t is the time, u the velocity field, P the pressure, 
v the viscosity and F the external force, acting at large 



scales. Turbulence is defined as a statistically stationary 
regime of the Navier-Stokes equations where an external 
forcing balances the viscous energy losses and induces a 
constant energy rate: du 2 jdt = cte. This means that the 
product F-u is scale independent, so that the force scales 
like 1/u. We note that eq. ^ exhibits common features 
with eq. 0. The l.h.s. of ([9| is invariant under the 
family of scale dilations <S/j(A) with arbitrary similarity 
exponent h and scale factor A [TU] : 

S h (X): £ — > X£, u^X h u, t -> A 1_h t, (10) 

where I is the space scale. This scale symmetry property 
is broken by both the forcing and the dissipation, which 
considered separately impose a dimensional scaling expo- 
nent: for the forcing, it is obtained by balancing F and 
du/dt, which given the transformation properties of F, 
yields u 3 ~ £, i.e. hp — 1/3. For the dissipation, balanc- 
ing of du/dt with v~S/ 2 u gives u ~ £ , i.e. hp, = — 1. 

Experimentally, none of these dimensional scalings is 
observed. Let us define the longitudinal velocity incre- 
ments at scale £ by 5ug — £■ [u(x+£)—u(x)]/£ and the ve- 
locity structure function of order n, S n (£) —< (Sui) n >, 
where <> is a statistical average. In a range of scales 
small compared with the forcing scale, and large com- 
pared with the scales at which dissipation becomes non- 
negligible (the so called inertial range), the structure 
functions are observed to follow approximate power law 
behavior: S n (£) oc £** n . The dimensional value im- 
posed by the forcing would be £n = nhp = n/3 while 
the dimensional value imposed by dissipation would be 
Cn = nfiD = —n. Experimentally, none of these dimen- 
sional scaling is observed. The origin of this anoma- 
lous scaling can be understood within the mechanism ex- 
plored previously: following the Navier-Stokes equation, 
the structure functions obey an equation of the type: 

I n S n = F n + D n , (11) 

where /„ is a scale invariant differential operator acting 
in the scale space, and F n and D n contributions due to 
forcing and dissipation, imposing dimensional exponents 
n/3 and —n. The scaling exponents in the inertial range 
then come from the scaling properties of the zero-mode of 
the operator I n . Note that until now, the explicit shape 
of the operator /„ has never been computed, even per- 
turbatively. It is quite possible that it is a non-diagonal 
operator of the structure functions, like in the Burgers 
equation, in which case the explicit computation of the 
zero- modes would be very difficult. Note also that for the 
Burgers equation, the scaling comes from the balance of 
non-diagonal contributions with terms corresponding to 
D n , which shows that the anomalous scaling described in 
the previous Section does not necessarily apply for any 
type of equation. In this respect, our model equation 
is a little bit ambiguous, and its application to turbulence 
not so straightforward as it may seem. We however adopt 
it as a working hypothesis, to explore its consequences. 
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Because of the presence of two scale-breaking symme- 
try mechanisms, turbulence offers a new possibility ab- 
sent in the model equation ([!]) : transition to new anoma- 
lous exponents at small scales, obtained by balancing of 
F and D. Because of constant energy transfers, F scales 
like 1/u, while D scales like u/£ 2 , so %t — 1. Indeed, ex- 
perimentally, one observes a transition towards S n oc £ n 
at small scales, the so-called "regular behavior" . The ex- 
istence of this transition imposes a scale dependence for 
the scaling exponents, and thus, non trivial scale depen- 
dence for the shape of the structure functions. Of course, 
this shape depends both on the shape of the operator /„, 
and on F n and D n . Given the symmetry properties of 
I n , it is however interesting to investigate whether the 
scale symmetry does impose possible generic shape for 
the structure functions, even within the inertial range, in 
the same way as it does constrain the n dependence of 
the scaling exponents. This is the purpose of the next 
Section. 



III. SCALE COVARIANT STRUCTURE 
FUNCTIONS IN TURBULENCE 

In the present Section, we assume that the velocity 
structure functions in turbulence obey a differential equa- 
tion which can be put under the form (111, where /„ 
is a scale invariant operator. Furthermore, we assume 
that the solution in the inertial range is given by the ho- 
mogeneous part of the equation (F n = D n — 0), and 
that there is a transition towards the "regular" solution 
S n (£) = K n £ n at a scale £ — r) n) where both K n and ij n 
are entirely determined by F n and D n . Our inertial range 
solution therefore extends up to rj n . Our aim is to derive 
a generic shape for the structure functions, satisfying all 
these requirements. From now on, we focus on stationary 
situations, thereby neglecting any time dependence. 



discrete slicing of the space scale I < L, under the form: 



A. On scale symmetry 

A scale invariant system is according to our definition, 
invariant by a transformation Sh(X) given in Note 
that in log- variables A = hxcf), T = ln£, this transfor- 
mation amounts to a translation invariance in A and T. 
This symmetry is rather intuitive, and connected with 
the possibility to multiply any characteristic quantity by 
an arbitrary constant, i.e. to perform arbitrary changes 
of units in the system [7j. This corresponds to a global 
scale invariance. In fact, as discussed by Pocheau [T2"] . 
scale symmetry encompasses another notion, correspond- 
ing to symmetry with respect to local scale dilations. 
This symmetry is connected with the invariance of the 
system with respect to changes of resolution, or equiva- 
lently to the possibility of performing arbitrary changes 
of rationalized (power law) unit systems [7J- To under- 
stand this local symmetry, it is convenient to consider a 



A 1 



L = T N £n 



(12) 



Here, T (or equivalently N) and A characterize the reso- 
lution of the observations. By local scale symmetry, we 
require that the system be invariant by changes of reso- 
lution, i.e. by changes of N. Changes in N are achieved 
by dilations of T: T — > T a , which transform i into ia. 
It is also a transformation which transform £i/£j into 
(£i/£j) a , and 4>i/<f>j into ((j>i/(j)j) a . A local scale dilation 
can therefore be seen as a transformation: 



£/£ -> (£/£ ) a , <t>/<j>o -> (0M>) a , 



(13) 



for arbitrary a, £q and 4>q. It can be seen as a change of 
rationalized unit system because it amounts to change £q 
and O by *o(W 1_a and MMfa) 1 "*- 

Summarizing, we say that a system is scale invariant, 
if it is invariant under both ^ and ( 13 ), i.e. if the equa- 
tions governing its evolution are covariant (keep the same 
shape) under these two transformations. 



B. Generic equation for structure function 

Assuming that a turbulent flow is scale invariant, we 
must then write the equations followed by a structure 
function S n {£) in a covariant way. For convenience, we 
introduce the log variables X n = In S n and T = In I 
[13] . The shape of the structure function is determined 
by the function X n (T), given some boundary (initial) 
conditions. We must write the differential equation fol- 
lowed by X n . Because of the global symmetry which 
amounts to a translation invariance in X n and T, this 
differential equation can only include derivatives of X n 
with respect to T . In the spirit of the amplitude equation 
theory, we write the generic equation as an expansion in 
power of the amplitude: 



dX n 
' dT 
dX n 
dT 



cPX„ 



d 3 X„ 



dT 2 n dT 3 
, dX„ d 2 X 

f fir, 



dT dT 2 



+ o{ 

0{{ 



d A X n 
dT A 

dX n . . 



dT 



)(14) 



where Xn,a n , b n , ■■■h n are some T independent constants 
characterizing the system and are to be chosen in order to 
respect the local scale symmetry. Applying the transfor- 
mation ( 13 ) to (141, we see that this equation is covariant 
only if the constants satisfy the following transformation 
rule: 



9n 



ah n ; . 



(15) 
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This suggests to introduce two characteristic constant 
scales in the system, £q and £\ and to write: 



Xn 
a n 
b„ 

fJn 

h„ 



Xn, 
Q"n. i 



= 6„ln(4/4), 
= c„(ln(V4)) 2 , 
= 9n, 



(16) 



where the constants £ n , ... are invariant under the local 
scale transformation, and only depend on the structure 
function (i.e. on n) and on the choice of £q and £\ . From 
now on, we fix the scale origin at the integral scale £ = L 
(representing the largest correlated motions). With this 
choice, scale invariance symmetry holds for 



ln(i?*) < T < 0, 



(17) 



where i?* = L/rj is a dimensionless number apparented 
to a Reynolds number and r\ is. the Kolmogorov scale 
(representing dissipative scales of turbulence). 



The case of turbulence 



Equation (14), with constants following (16) is a 



generic equation in a scale invariant system, satisfying 
both global and local symmetry. Its solution depends on 
the values of the constants, which must be computed fol- 
lowing some systematic procedure. However, interesting 
features of the solution may already be found in the lin- 
ear approximation, where all terms of degree A 2 or higher 
are neglected |14j . This approximation may be justified 
by considering a range of scales where the amplitude A 
is small enough. We study here only the simplest case 
relevant to turbulence. 



1. The linear case 

The highest relevant derivative is determined by the 
number of boundary conditions. The existence of a 
largest scale in the system (e.g. the injection scale) im- 
poses already at least one boundary condition. Further- 
more, the existence of the transition between the scale 
invariant solution and the regular solution at £ = r\ n im- 
poses at least a second boundary condition. Note that 
this is a peculiarity of turbulence. In many scale invariant 
systems, only one boundary condition is needed, corre- 
sponding to an infrared or ultraviolet cut-off in the scale 
space. 

We adopt L, the integral scale scale, and r/ n to nor- 
malize the constants, and introduce the pseudo- Reynolds 
number 



Dropping tildes for convenience, we therefore write the 
simplest differential equation relevant to structure func- 
tions in turbulence as: 



Xr, 



d\n{S n ) 
1 dT 



b n Itl(R T: 



d 2 HS n 
dT 2 



(18) 



where \n, o-n and b n are T independent constant 
depending on the order of the structure function. Note 
that they also are independent on the Reynolds number, 
since solutions with different Reynolds number R\ and 
i?2 can be related by a local scale transformation ( 13 ) 
with a = ln(i?!)/ ln(R 2 ). The solution is: 



Xn 



-T 2 + a n T + p n , a n = 0, 



HSn) - 2b n \n(Rn, 

ln(S n ) = ^T + a n e-^/^ R ^+(3 n , a n ^ffl) 
a n 

The constant a n and (3 n are fixed by the boundary con- 
ditions. Without loss of generality, we can fix the scale 
origin at £ — L, and choose the normalization for S n such 
that: 



ln(5„)(T = 0) 
d\n(S n ) 

jrp |T=-ln(i?.„) 



0. 



(20) 



The second boundary condition guarantees the transition 
towards the regular solution. These boundary conditions 
fix the constants as: 



a n = n+ /3 n = 0, a n = 0, 

On 



The solutions ( 19 ) with constant given by ( |21[ ) both fol- 
low the behavior schematized in Fig. [l] a power law 
behavior around £ = L, followed by a transition regime 
towards £ — rj n , where matching with the regular solution 
is obtained. 

The second solution, with a n ^ (generic case) is the 
most interesting, because it exhibits a feature reminiscent 
of what is observed in turbulence. This is the subject of 
the remaining Sections. 



2. Large Reynolds number limit 

In the limit where the pseudo- Reynolds number is very 
large, T/\n(R n ) 1 and the argument of the expo- 
nential can be expanded. The generic solution then ap- 
proaches the self-similar behavior: 



Xn 



(1- 



,-a n /b n 



)+ne- a " /K ))\n(£/L). (22) 



ln(i?„) = \n(L/j] r , 



This behavior is also valid in the neighborhood £ = L. 
This kind of asymptotic behavior has been observed in a 
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FIG. 1: The behavior of a structure function versus the scale 
when a n — (dotted line) and a n / (continuous line). 



tunnel wind experiment by Castaing et al |15] : when the 
Reynolds number is increased, the shape of the structure 
function goes from exponential power law, like in the first 
part of our generic solution (19d), to self-similar power 
law behavior. Even in the moderate Reynolds number 
regime, Castaing and his collaborators observe that the 
exponential power law regime is superseeded at larger 
scales by a pure power law regime. Such kind of behav- 
ior is indeed obtained at I > L in the generic solution, 
provided a n /b n is negative. This exponential power law 
as well as the variation of the exponent with (1/lni?) 
has also been predicted by Castaing [TH] using a scale in- 
variant Lagrangian formalism and requirement of finite 
dissipation. Within the present framework the result ap- 
pears to stem both from scale invariance and finite size 
effects (see discussion). 



3. Interpretation of the constants 

Up to this stage, we have not tried to characterize the 
various constants \n, o. n and b n appearing in the solu- 
tion. They should in principle be directly computed from 
the Navier-Stokes equations. Given the properties of the 
generic solution, we can however propose an interpreta- 
tion of the constants which provides some constraints on 
their shape. Consider the asymptotic behavior of the 
solution (22). It defines some "inertial range" scaling 
exponents: 



Cn 



= ne- a " /b " 



a n 



(23) 



They are made from two contributions: one stemming 
from the solution of the equation dA/dT = Xnl a n\ 
another one coming from the exponential power term, 
present only when b n is different from zero, i.e. when the 
system is of second order. As discussed in Section 3, the 
presence of this term is directly linked with the existence 
of a small scale cut-off, where the solution must bifurcate 



towards the regular solution. It can therefore be inter- 
preted as a finite-size effect. On the other hand, the con- 
tribution at b n = characterizes a perfect scale invariant 
solution, extending from to infinity in the scale space. 
This separation is reminiscent of the model of Dubrullc 
and Graner [7 , in which generic scaling exponents can 
be written: 



(24) 



Here, S( n is a contribution stemming only from scale 
symmetry considerations. It is characterized by its large 
n limit C = lim„_i. 00 S(^ n which can be interpreted as the 
codimension of the most intermittent structures. The 
term proportional to Aoo is the contribution due to the 
scaling properties of the maximum value of ln(S n ), which 
is defined only for finite-size systems. The factorization 



( 24 1 with a linear dependence of the contribution due to 



finite-size effect can also be justified from dynamical con- 
siderations [T7j. T his suggests to interpret the constants 



appearing in (231 as: 



-*oo : 



e -a n /b n = ^ 

Xn _ S( n 

cu 



1- A n 



(25) 



With this interpretation, finite size effects modifies the 
scale invariant solution in two ways: by the introduction 
of the linear term uAk, and by the modification of the 
codimension of the most intermittent structure: 



C s 



c 



finite 



(26) 



where C sca i e is the codimension in the scale invariant 
case, and C finite is the codimension in the finite size case. 
We stress that this finite size effect is Reynolds number 
independent (it only depends on the value of the constant 
b n ). It may however depend on other external parameter, 
such as the dimension (see discussion of Section 2.1). 
Using eq. (21), interpretation (26) completely deter- 



mines the ratio Xnl a n and a n /b n as a function of Aoo 
and the function <5£„, which can be severly constrained 
using only scale symmetry requirements [7]. 



4- Example: log-Poisson case 

As an example, let us considered the case where S( n 
corresponds to a log-Poisson statistics jT5H2"T] (one of the 
three possible cases derived in [7]): 



5c„ = c(i-/n 



(27) 



where C and /3 are two constants. If we adopt the val- 
ues C 'finite = 2, /3 3 = 2/3 and A m = 1/9 as advocated 
by She and Leveque |18j . for high Reynolds number tur- 
bulence, we can numerically determine the values of all 
parameters and thus, the scale dependence of the struc- 
ture functions, providing the pseudo-Reynolds number 
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n 




an/b n 


2 


0.53 


2.2 


3 


0.75 


2.2 


4 


0.95 


2.2 


5 


1.1 


2.2 


6 


1.25 


2.2 



TABLE I: Values of the parameters appearing in the structure 
functions for the log-Poisson model. 



R n are given. Table [I] summarizes the values obtained 
for n — 2 to 6. The corresponding structure functions, 
for R n = 10 are given in Fig. [2] 
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FIG. 2: The generic structure functions S2 to Sa in the log- 
Poisson case. 



5. Comparison with Batchelor's parametrization 

Using only symmetry arguments, we derived the shape 
of the structure functions in the linear approxima- 
tion. This shape can be compared with the Batche- 
lor's parametrization, which can be derived from the Kol- 
mogorov four-fifth law using matched asymptotics. This 
parametrization gives: 



1 + B(e/ V y< 



(2-C2V2' 



(28) 



where A and B are constants depending on the skewness 
of the velocity derivatives (B « 7.2 x 1(T 3 [22]). The best 
fit between Batchelor's parametrization and the generic 
solution is obtained for ry 2 = 6.35?7 and is shown in Fig. 

m 
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FIG. 3: Comparison between the Batchelor parametrization 
(dotted line) and the generic structure function h^S^) in the 
log-Poisson case (continuous line). 



6. General scaling 

The structure functions computed in the previous Sec- 
tion display an interesting property, referred to in tur- 
bulence as General Scaling (GS) [24]. We introduce the 
maximal event function Soc, defined as: 



Scv = lim 



S. 



n+l 



(29) 



In a bounded system, this function is always defined and 
represents the event of maximal intensity. In turbulence, 
it represents the largest velocity differences. The reduced 
structure functions Sn/S'!^ then obey a remarkable fac- 
torization property, 



In 



Sn \ S3 

$00 J V ^oo 



(30) 



This is shown in Fig. [4^. This factorization property, ex- 
tending throughout the whole scale interval, can be seen 
as a generalization of the factorization occurring in the 
inertial range, where all structure functions are propor- 
tional to T. 

This result can be explained by our choice of n in- 
dependent pseudo-Reynolds numbers R n = R. In such 
case, from (19) and p5|, we obtain: 



111(500) = A c 



ln(R) 



A T/ln(i?) _ 1 



(31) 



The reduced structure functions S n / S 1 ^ are then simply 
given by: 



In 



On 



1 



T 



ln(fl) 



(aW)_d 



(32) 

and are proportional to each other within the whole scale 
interval. It is not clear whether the present analysis pro- 
vides an explanation to the GS observed in turbulence. It 
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solution occurs at the same scale rj n = f] (see [23] for 
a discussion about n dependent cut-offs). If the depen- 
dence of ln(i? n ) with n is weak, as seems to be the case 
in turbulence, then the present analysis is still valid in 
turbulence and explains the phenomenon of GS within 
the context of scale symmetry. 



7. Extended Self-Similarity 



The factorization property ( 30 ) can be strengthened 
if, in addition, ln(5oo) is also proportional to any 
ln^n/S^). In that case [24], the logarithms of any struc- 
ture functions are proportional to each other. In other 
words, when one structure function is plotted against 
another one (e.g. the third one), one observes a well 
defined scaling regime, even in the range of scale where 
the function is not self-similar. This property was called 
Extended Self-Similarity [53]. The condition to observe 
this is Aqo — > 0. In practise, if this condition is satis- 
fied approximately, one can expect to observe ESS in the 
system. In our case, Aoo = 1/9 is small, so we should ob- 
serve ESS. This is illustrated in Fig. \4)p. Note that ESS 
property means that the relative exponent £* defined as 



C* 



dln(S n ) 



dln(S 3 ) ' 

is much better defined than the true exponent: 



dhi(S n ) 
d\n(£) 



(33) 



(34) 



To illustrate this point, we have computed these scaling 
exponents for n — 6 (Fig. |4j;). It can be seen that ^6 
decreases steadily from £ 6 = 6 to £ 6 = 1.78, while Q 
displays much weaker variations over the whole interval, 
from (g = 2 to Q w 1.8 (the "intermittent value"). 



8. Finite size effects vs asymptotic K41 solution? 



FIG. 4: Top: (a) General Scaling: all the reduced structure 
functions are proportional to each other (here the third one). 
Middle: (b) Extended Self Similarity: when the 6 th order 
structure function is plotted against the 3 r one, one observes 
a better defined scaling regime (compare with Fig. 2). De- 
spite the visual impression, the Extended Self Similarity is 
only approximate (see Section 3). Bottom: (c) Local scal- 
ing exponent fg and relative scaling exponent C6/C3 versus 
the scale. The local scaling exponent decreases constantly 
with the scale, proving the absence of exact self-similarity. 
By contrast, the relative scaling exponent varies much less, 
which enables a better determination of this exponent (ESS 
property) . 



is not known in turbulence whether the pseudo-reynolds 
numbers R n are actually n independent, i.e. if the tran- 
sition from inertial range solution towards the regular 



In the previous Sections, we have used an interpreta- 
tion of the constants to compute explicitly structure func- 
tions. This interpretation was dictated by our choice to 
introduce explicitly the dissipative range in the boundary 
condition (via the matching to the "regular solution"). 
In such interpretation, asymptotic (high Reynolds num- 
ber) scaling exponents take the shape ( 23 1 , with a lin- 
ear part coming explicitly from finite size effects. In ab- 
sence of finite size effects, they take the simple shape 
C« = C sca i e (l — f3 n ) predicted by Dubrulle and Graner 
[7]. The Kolmogorov solution £ n = n/3 is never reached, 
unless the coefficient (3 appearing in ([19]d) depends on 
other external parameters such as the dimension. The 
Kolmogorov solution could then appear as the infinite 
dimension limit (even in presence of finite size effects), 
and j3 would be a parameter characterizing the dimen- 
sion of the system. 
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Other interpretations are however possible, within the 
same symmetry arguments. For exemple, one could re- 
strict the solutions to an "inertial range" of scale defined 
by imposing (3 = 1, i.e. 0:3 = in (27). Asymptotic 



Kolmogorov solution could then be obtained with the 
following choice of constants (still compatible with the 
theory of Dubrulle-Graner): 



Xr 

Or, 



-Or. 



b n ln(R n ) 

p -a n /b n 



S_Cn 

n 8( n 
3 6( 3 - 

3Aoo. 



(35) 



In such case, the derivative dln(S n ) / dT reaches the value 
(24) at T = ln(i?„), which appears as another bound- 



ary condition. In absence of finite size effects, the scal- 
ing exponents are £„ = (1 — /3 Tl )/(l — /? 3 )- The Kol- 
mogorov solution Cn/3 is obtained as the asymptotic so- 
lution (Reynolds tends to infinity) in presence of finite 
size effects. 

We were not able to discriminate between these two 
types of interpretations. Obviously, one is rather valid 
in the vicinity of the dissipative range. It could then be 
seen as a refinement of the arguments by Frisch and Ver- 
gassola [25] obtained within the multifractal model (see 
|15j for a discussion within the log-similarity hypothe- 
sis) . The second is valid in the "inertial range" , defined 
using £3 = 1 (inertial range log-similarity hypothesis, see 
|15j). They however lead to distinct prediction about the 
possibility to approach Kolmogorov solution. The most 
recent experimental results seem to indicate that the ob- 
served scaling exponents are almost independent of the 
Reynolds number This would favor the first inter- 
pretation, and explain our choice in the present paper. 
Obviously, it would be interesting to consider further the 
variation of scaling exponents with the Reynolds number, 
and, possibly, to reconsider the second interpretation. 

It should be stressed however that in any case, the 
"exponential power-law" appears as a prediction of sym- 
metry arguments, independent of any boundary consid- 
eration, i.e. of these kind of interpretations. 



scale saturation of the structure functions, which is ab- 
sent in the linear model, could be obtained by taking into 
account non-linear terms in the differential equation, i.e. 
by a slightly more complicated model. 

Our results illustrate the essential influence of scale 
symmetry on structure functions in turbulence, and pro- 
vides further support to the scenario of " scale invariant" 
anomalous scaling discussed in Section 2. We note that 
finite size effects (ultraviolet cut-off) generically lead to 
non-power law behavior of the structure function, but 
rather to exponential power-law behavior, exp(£ a ), where 
a is real, proportional to the inverse Reynolds number. 
Such dependence, connected with the requirement of co- 
variance by resolution, was also inferred by Barenblatt 
and Goldenfeld using the principle of " Reynolds number 
covariance" j^Hj- We see here that such principle (left 
unjustified by Barenblatt and Goldenfeld) is directly con- 
nected with scale symmetry. Finally, we note that there 
is a possibility to get complex exponential power-law be- 
haviors (a complex) if we allow a differential equation 
of higher order, or if we allow the presence of terms di- 
rectly proportional to A in (14). The first possibility 



could be justified if more than two boundary conditions 
are necessary to specify the solution. The second pos- 
sibility still requires only two boundary conditions, but 
implies a breaking of global scale invariance. This re- 
quires the existence of a privileged scale into the sys- 
tem, and occurs for example in a system subject only to 
discrete scale invariance |27j . Complex exponents gives 
rise to log-periodic oscillations at large Reynolds number 
(when the solution goes from exponential power law to 
power law), which may have been detected in a variety of 
physical systems [27l EE] • It would be interesting to see 
whether they can also arise in certain turbulent flows. 
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IV. DISCUSSION 



Using only symmetry considerations, we were able to 
build generic structure functions reproducing many fea- 
tures observed in actual structure functions in turbu- 
lence: transition from exponential power shape to a 
power shape with increasing Reynolds number, extended- 
self-similarity, regular matching with the regular solution 
at small scale. We do not claim that the generic struc- 
ture functions considered in the present paper actually 
fit exactly the structure functions determined experimen- 
tally, because we intentionally considered only the sim- 
plest case relevant to turbulence, where the equation is 
linear. It is clear for example that the observed large 
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